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One of the puzzling aspects of high temperature 
superconductors is the prevalence of magnetism 
in the normal state and the persistence of super- 
conductivity in very high magnetic fields. Gen- 
erally, superconductivity and magnetism are not 
compatible. But recent neutron scattering results 
indicate that antiferromagnetism can appear 
deep in the superconducting state in an applied 
magnetic field. Magnetic fields penetrate a su- 
perconductor in the form of quantized flux lines 
each one representing a vortex of supercurrents. 
Superconductivity is suppressed in the core of the 
vortex and it has been suggested that antiferro- 
magnetism might develop there [Q. To address 
this question it is important to perform electronic 
structural studies with spatial resolution. Here 
we report on implementation of a high field NMR 
imaging experiment [^|^ that allows spatial res- 
olution of the electronic behavior both inside and 
outside the vortex cores. Outside we find strong 
antiferromagnetic fluctuations, and localized in- 
side there are electronic states rather different 
from those found in conventional superconduc- 
tors. 

Supercurrents that form a vortex are expected to decay 
inversely with distance from the vortex core giving rise to 
a spatial distribution of internal magnetic fields, P{IIi„t), 
shown in Fig. |^. The color coding identifies spatial po- 
sitions in the lattice of vortices with their corresponding 
location in the field distribution for a vortex lattice. This 
field distribution will be the same as the NMR spectrum 
for any nucleus with small intrinsic broadening as is the 
case for ^^O in YBCO. Its narrow spectrum in the nor- 
mal state (blue trace) is shown in the figure for two of 
the four quadrupolar satellite resonances and compared 
with the vortex broadened spectrum (red trace) in the su- 
perconducting state 1^,0. We isolate the —1/2 ^ —3/2 
transition giving the black trace, corresponding to the 
expected field distribution shown in the inset. Conse- 
quently, our NMR experiment allows a spatially resolved 
study of the vortex structure. During the 
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FIG. 1. Planar spectra versus the internal mag- 
netic field. Hint, showing broadening from vortices. 
Hint = '^/^'^'y — Ho where ^^7 is the gyromagnetic ra- 
tio for oxygen and the spectrometer frequency, u>, is set 
to resonance at the pealc frequency for the -1/2 ^ -3/2 
transition in the applied field Ho=13 T at 11 K. For clar- 
ity only two quadupolar satellite resonances are shown. 
In the normal state at 100 K (blue) there is a sharp edge 
which broadens at low temperature (red). The spectra 
are shifted such that internal fields are measured relative 
to the peaic. The black trace represents the -1/2 <-> -3/2 
part of the spectrum obtained by subtraction of the 
-3/2 <-» -5/2 transition. The latter was obtained from the 
convolution of the normal state spectrum with a broaden- 
ing function measured independently on the 3/2 ^ 1/2 
transition (not shown here). Inset: We have calculated 
the internal magnetic field profile and the correspond- 
ing field distribution function for a 80° vortex lattice at 
Ho = l3 T, with a penetration depth of A = 1500 A. The 
minimum field at the center of the vortex lattice unit 
cell (m), the maximum field at the vortex cores (M), and 
the saddle-point field midway between vortices (s), are 
shown. 

NMR experiment the oxygen nuclei are perturbed and 
then relax back to thermal equilibrium at the spin- lattice 
relaxation rate, which is determined by the density of 
electronic states near the Fermi energy. We have made 
spin-lattice relaxation rate measurements, Fig. ^, that 
are also spatially resolved and we interpret our results 
in terms of the electronic excitation spectrum both in- 
side and outside the vortex cores. Within the core we 
expect NMR to be sensitive to localized electronic states 
such as have been observed in conventional superconduc- 
tors and recently in high temperature superconductors 
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1^,0 . High magnetic fields allow us to use the Zeeman 
effect to study such electronic excitations. Varying the 
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FIG. 2. Planar ^^O spin-lattice relaxation rate across 
the vortex lattice NMR spectrum at 11 K. The spectrum 
is the -1/2 « -3/2 transition at 37 T. The shaded re- 
gion is the estimated area occupied by the vortex cores 
at 37 T. We used a superconducting magnet, 13 T; a Bit- 
ter-style electromagnet, 23 T; and a hybrid of the two 
technologies, 37 T. We obtained precise spectra using a 
field sweep technique and we measured the relaxation rate 
using progressive saturation p^ . 

field up to 37 T, we sweep through a significant range 
of corresponding Zeeman energies, ~ ±2.5 meV, larger 
than the thermal energy. 

Our aligned powder sample of YBa2Cu307_a is near- 
optimally doped and ^ 60% ^''0-enriched. The crystal 
c-axis was aligned with the direction of the applied mag- 
netic field. Low-field magnetization data show a sharp 
transition at Tc(0) = 92.5 K. High magnetic fields are es- 
sential to increase sensitivity to vortex cores by increasing 
their fraction of the NMR spectrum. In a field of 37 T 
vortices are ~ 80 A apart as compared with their core 
radius given by the size of the electron-pair wavefunc- 
tion, ^ « 16 A, and so the cores occupy ^ 15% of both 
the total volume of the sample and the NMR spectrum. 
We have shaded the corresponding region in Fig. ||. The 
high field also suppresses spin diffusion and vortex vibra- 
tions which complicate interpretation of the spin lattice 
relaxation at low field 0. 

For high-temperature superconductors the gap in the 
excitation spectrum, Ak, is believed to have d-wave sym- 
metry, implying that there are 4 nodes on the Fermi sur- 
face where the gap vanishes. At low temperature the 
electronic excitations, called quasiparticles, come exclu- 



sively from the nodal regions shown in Fig. They can 
be probed with NMR spin lattice relaxation, with the 
rate given by, 
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FIG. 3. Sketch of the Fermi surface for YBCO and 
the three processes that contribute to the NMR rate at 
low temperatures. The three processes, Fi, F2, and F3, 
correspond to the choice of nodes, points where super- 
conductivity is fully suppressed, (black dots) for the ini- 
tial and final quasiparticle states. The spin lattice re- 
laxation rate is given by, T^^ ~ {Ni{e)N f{e)) , where 
No. = X^k ~ i?k,a) are the initial (a = i) and final 
(q = /) density of states (DOS) respectively. The sum 
over wavevector k is restricted to the regions around the 
nodes, k^, and the product of the DOS is integrated over 
energy e in the range ksT near the Fermi energy, Ep. The 
excitation energies i?k,a for quasiparticles moving in the 
superfiow field with momentum pa axe Doppler shifted 
by an amount Da = v_F,a • Ps, where the Fermi velocity 
at node k^ is given in terms of the normal state excita- 
tion energy tk by, vf.q ~ 9ek/9k|k=ka. Note that for 
process F3 the relation Di — —Dj holds, in contrast to, 
for example, process Fi where Di — Df. In sufficiently 
high magnetic field the Zeeman energy, becomes impor- 
tant and enters initial and final stat es with d ifferent signs 
{Z,,f = ±Z), leading to, Ek.a = + ^k - + Da,. 
Since the DOS depends linearly on energy near the nodes, 
the NMR rate is proportional to the product given by 
Eq. I 

T,"'^{\e-Z + D,\\e + Z + Df\), (1) 

where e is of the order of fc^T; Z = —^je^iHo, is the 
Zeeman energy; Da — ^p.a ■ Ps is the Doppler shift, and 
a = i, f corresponds to initial and final electronic states 
near the nodes. Depending on the relative magnitude of 
ksT, Z, and D the dependence of the relaxation rate 
on temperature and field will show different behavior for 
each of three processes shown in Fig. ^ We can vary 
the first two terms, ksT, Z, by changing temperature 
and magnetic field. The Doppler term is proportional to 
the vortex currents, and varies with internal field. Hint, 
across the NMR spectrum. This term is a microscopic 
analogue of the classical frequency shift of a wave in a 
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moving reference frame, the latter being associated with 
the superconducting condensate. 

At the peak of the NMR spectrum. Hint = in Fig. ||, 
we observe that the rate increases as a function of the 
external field, Hg. This dependence comes from the Zee- 
man interaction. The Doppler term is small in this re- 
gion since the supercurrents cancel at the saddle-point. 
As we increase the applied field the Zeeman interaction 
dominates the energy spectrum near the nodes, i.e. Z, 
equivalent to 25 K at 37 T, is larger than ksT at 11 K, 
and from Eq. Ill the NMR rate rises quadratically with 
field as (1^2 -P|). 

To the right of the peak in the spectrum toward the 
vortex core. Fig. ||, we find that the rate increases by a 
factor of ~ 4—5. Here the superflow momentum increases 
sufficiently that the Doppler term in Eq. ^ becomes large 
compared to both Z and fc^T and it dominates the quasi- 
particle energy. Our observations show that T^^, and 
consequently the density of electronic states (DOS), de- 
pends on distance, r, from the core. In particular, we 
expect the rate to increase as T^^ ^ (ps) ^ r^^. The 
enhancement of T^^ depends only on r and should not 
be dependent on the applied magnetic field other than 
the Zeeman interaction. This is not to be confused 
with Volovik's prediction ||ll| that the average DOS, 
(iV(0)) y^Ho, which corresponds to (Tf ^) HglnHo 
where the average is over the entire spectrum. For con- 
ventional superconductors without nodes, the Doppler 
term has a negligible effect. Our observation of the 
Doppler shift extending to the saddle-point region shows 
that, even in high magnetic fields, the superconducting 
gap has well defined nodes. 

Outside the vortex core we find that the effect of ap- 
plied magnetic field from 13 to 37 T is simply to displace 
the relaxation rate data on a vertical scale in Fig. ^ This 
is possible if the Doppler term in Eq. |l| changes sign for 
initial and final states, just like the Zeeman term, and 
can only occur if those quasiparticle states are not from 
the same node. Then, at low temperatures, the rate is 
~ + Z"^ and increases quadratically with field. Oth- 
erwise, initial and final states come from the same node 
and the rate is ^ Z?^ — Z^ , having the wrong field depen- 
dence in the regime D > Z . We have calculated the rel- 
ative strengths of the three possible processes and their 
dependence on magnetic field. The calculation is per- 
formed using the known geometric form factors and the 
imaginary part of the electronic susceptibility which we 
extract by fitting the data. To account for our observed 
field dependence, the spin susceptibility must be strongly 
peaked near a wavevector corresponding to the F2 or F3 
process, at least a factor of 7 greater than that at zero 
wavevector. This indicates that strong antiferromagnetic 
correlations between quasiparticles coexist with super- 
conductivity at low temperatures. Antiferromagnetically 
correlated quasiparticles were identified by neutron scat- 
tering ||l| in LSCO, a related compound, and are now 
confirmed in YBCO from NMR. Lake et al. speculated 
that if there was antiferromagnetism in the vortex core it 



might polarize the intervening medium and account for 
their experiment. 

Near the vortex core, the NMR rate decreases after 
reaching a maximum in the applied field of 13 T as shown 
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FIG. 4. Spin-lattice relaxation rate of planar ^^O di- 
vided by temperature as a function of internal magnetic 
field. The shaded region is the area occupied by the vor- 
tex cores at 37 T. For a d-wave superconductor in the 
temperature range where ksT is much less than yehHo, 
{TiT)~^ is expected to be constant. We clearly observe 
this behavior near the peak in the NMR spectrum. 

in Fig. H, refiecting a drop-off in p^. At higher applied 
fields the rate increases significantly, well beyond its max- 
imum value at 13 T. This increase with field exceeds that 
measured outside the cores by a factor of four and is most 
pronounced at the lowest temperatures. Fig. ^. We at- 
tribute this substantial excess rate to vortex core states. 
When superconductivity is suppressed near surfaces or 
in the vortex core, localized electronic states, called An- 
dreev bound states, can appear P,p^. Their signature is 
a peak in the density of states at the Fermi energy, called 
a zero-bias anomaly in tunneling experiments. Simi- 
lar behavior might be expected in the vortex core of a 
high temperature superconductor. If we interpret our 
relaxation rate in the same way as for quasiparticles out- 
side the cores, following the discussion in the caption to 
Fig. H and Eq. |[ then both field and temperature depen- 
dence are opposite to what is required for the zero-bias 
anomaly. Consequently, instead of a peak in the density 
of states there is a mini-gap, ~ ±5 meV, much sharper 
than the variation of the DOS outside the vortex core. 
Previous reports of vortex core states have been based on 
charge tunnehng in YBCO [[s) and BSSCO Near 
Ep, ~ ±5.5 meV for YBCO, an added structure in the 
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DOS from the vortex core region separated by a mini- 
gap was observed. This surprising result is inconsistent 
with theoretical predictions Furthermore, tunnel- 

ing is extremely sensitive to the surface and depends on 
an understanding of the tunneling matrix elements p6| . 
In contrast NMR probes bulk material. And yet, both 
experiments can have a similar interpretation. 

NMR is directly sensitive to the magnetic character of 
electronic excitations. So an alternate explanation of our 
results for the vortex core region might involve antifer- 
romagnetism as predicted [2| from a S0(5) nonlinear a 
model. In this case field induced antiferromagnetic exci- 
tations could strongly enhance the NMR rate consistent 
with our observations. In future work, spatially resolved 
NMR relaxation experiments can be extended to under- 
stand the effects of doping on vortex core states and pos- 
sibly identify a connection with the spin-pseudogap found 
in the normal state. 
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